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Recent theory and experiment in crystals of molecular magnets suggest that 
fundamental tests of the decoherence mechanisms of macroscopic quantum 
phenomena may be feasible in these systems (which are almost ideal quantum 
spin glasses). We review the results, and suggest new experiments. 

q The attempt to push quantum effects to mesoscopic or macroscopic 

O ■ scales has led, in magnetic systems, to various interesting discoveries. The- 



oretical work on domain wall tunnelinga'Oa led to experiments claiming to 
^ \ see tunneling of single domain walls Jj Early theory on tunneling of large 

' spinsi' i stimulated e xperime nts on resonant tunneling in crystals of magnetic 

macromoleculesP&ollololl 2 ] and on "macroscopic quantum coherence"" in 
giant ferritin molecules. 

There is a theoretical complication - the environmental "spin bath" 
of nuclear and paramagnetic spins should destructively interfere with any 
tunneling. BEER In fact, except in ferritin, all experiments see incnheiepi 
tunneling relaxation at low T. A detailed theory of such relaxation, Jla>llH>E2'E!l 
applied to ensembles of tunneling nanomagnets, makes several "universal" 
predictions for the relaxation characteristics, some of which have now been 
verifiedSc 3 ! Here we discuss the recent developments, and suggest future 
experiments designed to bring out the fundamental role of nuclear spins in 
quantum decoherence at low T. 
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1. Kinetics of Quantum Relaxation 

Our problem is simple but general. Consider an array (ordered or oth- 
erwise) of identical 2-level systems (which here represent large spins; S S> 1, 
below the crossover temperature T c of confinement to their 2 lowest levels). 
Their Hamiltonian is" 

H = \ E v ^^ ] + E ^ ] + E {h}) + h nn , (i) 

3 jk 

where V-j describes dipolar interactions between nanomagnets {i,j}, A is 
the tunneling amplitude between states a = ±, V^r couples nanomagnets to 
nuclear spins {Ik}, and Hn n describes the intrinsic nuclear spin dynamics^ 
Note that without Vn, Eq. (|l|) describes the "quantum spin glass" problem;! 2 ^ 
the effect of Vm is to couple this to an environment. 

The entire system is described by normalised probability distributions 
for single molecules, pairs of molecules, etc.; thus P a (£, r; t) is the probability 
that a molecule at position r has polarisation a = ± and bias £ = e+ — e_ 
between states |±) = | ± S). The bias £ at r = fl sums the nuclear bias 
(of typical energy scale Eq), the dipolar bias J2j VyT^ (typical energy scale 
Ed), and any bias from external fields. This 1-particle distribution Pry 
is connected to higher distributions P( 2 )>-P(3)) e t c -i y i a a BBGKY- style 
hierarchy of kinetic equations.^ Since the only part of the environment 
with any dynamics at low T and low H (ie., able to mediate irreversible 
incoherent tunneling) is the nuclear bath, the kinetic equation for Pry is 
governed by a nuclear-mediated transition rate tZ 1 , calculated inJll When 
the number of nuclear spins co-flipping with the molecular spin is small, 

TnHO ~ T^e-WC , (2) 

where £ G parametrises the range of fluctuations in £, and t^ 1 ~ A 2 /£ Q ; 
typically A < £, < E , E D . 

Several unambiguous predictionsEH flow from the kinetic equation. Over 
a wide time range (Ed/^q > t/ro > £q/Ed) the magnetisation should re- 
lax from saturation according to a universal square-root law, i.e., M(t) = 
M [l - (t/TQ) 1 ' 2 }. For ellipsoidal samples r Q t$ = (E D T /£ )f(c) 
(where /(c) is an analytic function of the ellipsoid shape); for other shapes 
t q ~ r^"/£ -/V T (0), where JV T (0) = /d 3 rP T (f = 0,r,i = 0) is the initial "up" 
distribution at zero bias. Since Tq ~ Ed /A 2 , we have tq ~ E d /A 2 S,q for a 
generic shape (since iVj(O) ~ 0(1/ Ed))- In Figure 1 we see why the square 
root prediction is valid independently of the details of the hyperfine cou- 
plings, the sample shape, or the ratio Eq/Ed', it arises from the "Lorentzian 
hole", eaten in the initial iVt(£) by the decay, coming from the long-range 



Quantum Relaxation of Ensembles of Nanomagnets 
0.005 




CO 
CD 



-0.005 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

-300 -200 -100 100 200 300 

Energy 

Fig. 1. 

dipolar inter actions. Fig. 1 shows a Monte Carlo simulation for a cube of 
size (75) 3 sites, and depicts A T (£) when t = (so M = M ), when M/M has 
relaxed to 0.99; and the difference between the two. In a weak external field, 
we must replace iVj(O) in the expression for Tq by iVf(£ = g/JBSH); thus 
a field scan will find that Tq 1 oc iVj(£ = gusSH), ie., Tq (H) directly mea- 
sures iVf(£). Since iVj-(£) can be calculated numerically (knowing Eq,Ed, 
and the sample shape), this is a testable prediction. 



2. Universality and Crossovers 



How "universal" is the square root law? One potential problem, given 
that the ligand groups in such large molecules can be unstable, is rogue 
fast-relaxing "impurity" molecules. Suppose the impurity concentration is 
5 < 1. Defining M A (t) = (1 - S ) - x(t), and M B (t) = 5 - y(t), for the 
bulk and impurity magnetisations, normalising Mq = 1, and assuming a ratio 
A > 1 between impurity and bulk relaxation rates, one gets a short-time 
behaviour 



y(t) = <J [1 - exp (-Ax(t)/(1 + <J ))] 



(3) 



t/rA = x 2 (t) + 28 



x(t) H — I exp 



-A 



(4) 
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One easily verifies that for arbitrary A5q, any appreciable impurity con- 
tribution to dM/dt destroys square-root relaxation (conversely, square-root 
relaxation only appears if the impurity contribution is negligible). 

The square root behaviour must also break down in the crossovers to 
long-time relaxation, and to thermal activation. After ~ 10 — 15% of the 
sample has relaxed, the Lorentzian hole will distort because of intermolecular 
correlations. These are not because of multi-flip processes (which have rate ~ 
A^/^qEd, and so are rare), but because factorisability of fails. However 
at long times other universal results emerge. First, if hyperfine spread in 
-ZVt(£) is strong (ie., when Eq > Ed), then at times > Ti(T), the T\ nuclear 
relaxation will sweep at any molecule over most of the range A£ of N(£); 
one then getsa M(t) ~ er t / Tl ^ T \ for t 3> T%. For weak hyperfine coupling 
{Eq <C Ed), intermolecular correlations give a complicated quantum spin- 
glass problem, with no "sweeping" mechanism to wipe out these correlations. 
Second, if we thermally recycle the sample, depolarising at high T until 
M(t)/Mo <C 1, the remaining relaxation will be exponential, at a rate0 
(with k ~ 1): 

t E d [1 + k\h(Ed/^ )} 

The quantum-thermal crossover is more subtle. The spin gap in the 
system is Auj = (E$ — £L(s_i)), between the lowest levels | ± S) and the 
next 2 levels | ± (S — 1)), and we expect that roughly kT c ~ Acj/27t (in fact 
Auj ~ 11.6-ftT (5. IK) for Mn-12 (Fe-8) respectively, whereas experimentally 
kT c ~ 2K (OAK). As noted above, if T T c , everything but the nuclear 
subsystem is frozen (dipolar flip-flop processes to levels ^(S 1 — 1)) occur with 
frequency ~ Q.di p e~^ u) f kT , where the coupling Sldip ~ 10 6 — l$Hz in these 
systems; at the lowest temperatures (70mK) used in the Fe-8 experiments,!!! 
the activation exponent was ~ e -70 !!). In this case we expect square root 
relaxation; and tq will only depend on temperature to the extent that Ti 
does (ie., very little). However as one approaches T c , molecular flip-flops 
are more frequent. The square root prediction is not affected- but tq can 
develop a thermally-activated component, because (a) now dipolar field fluc- 
tuations can also bring molecules to resonance, and (b) by exciting nuclear 
T\ transitions, they sweep the nuclear bias (this mechanism will only be im- 
portant for strong hyperfine interactions). A detailed theory of this is quite 
complicated.il 

Finally, in the very low-T limit 2 things will happen; the nuclear spins 
will polarise in the hyperfine field, for temperatures below the hyperfine 
coupling uj), and the nuclear T2 fluctuations will freeze out; and glass theory 
then predicts that a "dipolar gap" will open up in N a (^) if M Mq. 
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3. Experiments 

Some of the predictions (such as the square root law in the quantum 
regime) have recently been found, in experiments on Fe-8 and Mn-12 single 
crystals .El El However there is still no direct proof for the controlling role of 
nuclear spins in the quantum dynamics of nanomagnets; the observed \Jt/rQ 
relaxation clearly shows the role of dipole interactions, but the evidence for 
the nuclear role is indirect. How can it be seen directly? 

One obvious way is to do AC absorption experiments over a frequency 
range encompassing both A and £ Q . In Mn-12 and Fe-8, A is perhaps 10~ 10 — 
10 -8 -fT, some 8-10 orders of magnitude smaller than Ed or Eq (which is 
why it is so astonishing that one sees very low T relaxation at all); we 
expect £ G ~ W~ 5 K. This test is discussed by Rose and Stamp.c3 A much 
more intriguing possibility would be available if NMR experimentalists could 
find resonance lines for some of the spin-5/2 Mn nuclei in Mn-12 (the zero- 
field hyperfine coupling is probably ~ 600 MHz (450 MHz) on the spin-2 
(3/2) Mn electronic sites respectively). At this point a "David and Goliath" 
experiment becomes possible; one "tickles" the nuclear spins at the resonance 
frequency, in order to control the relaxation of the giant magnetic molecules 
(eg. retarding the relaxation by changing £ D ). Another, probably more 
difficult experiment would "freeze out" the nuclear dynamics by cooling to 
kT <C uj (ie, well below 30mK for Mn-12). The problem here is that of 
cooling the nuclei, since most of the heat transfer will involve the same 
T2 fluctuations that we want to freeze out! However there is a fundamental 
interest to doing this, since according to the theory these nuclear fluctuations 
are the last barrier to truly macroscopic coherent (MQC) behaviour in low-T 
systems (note that recent workE3 indicates that nuclear and paramagnetic 
spins will make it very hard to see MQC even in superconductors; spin baths 
tend to be far more damaging to MQC than oscillator bathsc3). This, as well 
as the resulting "Quantum Spin Glass" behaviour in the unpolarised limit, 
means that the low-T limit promises much to the intrepid experimentalist! 
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